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Abstract 

We show how to reahze a "protected" qubit by using a fully frustrated 
Josephson Junction ladder (JJL) with Mobius boundary conditions. Such 
a system has been recently studied within a twisted conformal field theory 
(CFT) approach [I112] and shown to develop the phenomenon of flux fraction- 
alization [3]. The relevance of a "closed" geometry has been fully exploited 
in relating the topological properties of the ground state of the system to 
the presence of half flux quanta and the emergence of a topological order has 
been predicted In this letter the stability and transformation properties 
of the ground states under adiabatic magnetic flux change are analyzed and 
the deep consequences on the realization of a solid state qubit, protected from 
decoherence, are presented. 
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Arrays of weakly coupled Josephson junctions provide an experimental realization of 
the two dimensional {2D) XY model. A Josephson junction ladder (JJL) is the sim- 
plest quasi-one dimensional version of an array in a magnetic field [5J; recently such a 
system has been the subject of many investigations because of its possibility to display 
different transitions as a function of the magnetic field, temperature, disorder, quantum 
fluctuations and dissipation. In a recent paper [5] we analyzed the phenomenon of frac- 
tionalization of the flux quantum in a fully frustrated JJL in order to investigate how 
the phenomenon of Cooper pair condensation could cope with properties of charge (flux) 
fractionalization, typical of a low dimensional system with a discrete Z2 symmetry. The 
role of such a symmetry was recognized to be crucial for demanding more general bound- 
ary conditions, of the Mobius type, at the end sites of the ladder [3]. The same feature 
was evidenced also in quantum Hall systems in the presence of impurities or defects [B] 
[7| [8J . Furthermore a Z2 symmetry is present in the fully frustrated XY (FFXY) model or 
equivalently, see Refs. [9] [ID], in two dimensional Josephson junction arrays (JJA) with 
half flux quantum threading each square cell and accounts for the degeneracy of the 
ground state. We noticed how it was possible to generate non trivial topologies, i.e. the 
torus, in the context of a CFT approach, which allowed us to construct a ground state 
wave function, whose center of charge could describe a coherent superposition of localized 
states sharing all the non trivial global properties of the order parameter. In particular 
for the FFXY model they were shown to be closely related to the presence of half flux 
quanta, also viewed as topological defects[3]. The emergence of topological order in fully 
frustrated JJLs with non trivial geometry has been predicted and fully exploited in Ref. 
[1] by means of CFT techniques. Such a concept was first introduced in order to describe 
the ground state of a quantum Hall fluid [llj but today it is of much more general interest 
|12] . Two features of topological order are very striking: fractionally charged quasiparti- 
cles and a ground state degeneracy depending on the topology of the underlying manifold, 
which is lifted by quasiparticles tunneling processes. In general a system is in a topolog- 
ical phase if its low-energy, long-distance effective field theory is a topological quantum 
field theory that is, if all of its physical correlation functions are topologically invariant 
up to corrections of the form e~ r at temperature T for some nonzero energy gap A. More 
recently superconductors have been proposed in which superconductivity arises from a 
topological mechanism rather than from a Ginzburg-Landau paradigm: the key feature is 
a mapping on an effective Chern- Simons gauge theory, which turns out to be exact in the 
case of JJA and frustrated JJA [13j. As we will stress in the following, topological order 
is crucial for the implementation of fully frustrated JJLs as "protected" qubits [H] [15] in 
solid state quantum computation realm. The idea in all such realizations is that the sys- 
tems involved (large and small size Josephson junction arrays of special geometry [16]) 
share the property that, in the classical limit for the local superconducting variables, the 
ground state is highly degenerate. The residual quantum processes within such a low en- 
ergy subspace lift the classical degeneracy in favor of macroscopic coherent superpositions 
of classical ground states [_16j. An example of such a system has been proposed, which 
consists of chains of rhombi frustrated by an half flux quantum [16] with the property 
that in the classical limit each rhombus has two degenerate states. The protected degen- 



2 



eracy in all such systems emerges as a natural property of the lattice Chern-Simons gauge 
theories which describe them [16]. In general, if a physical system has topological degrees 
of freedom that are insensitive to local perturbations (that is noise), then information 
contained in those degrees of freedom would be automatically protected against errors 
caused by local interactions with the environment |15j . 

The aim of this letter is to show how to realize a "protected" qubit in terms of a fully 
frustrated Josephson Junction ladder (JJL) with Mobius boundary conditions by fully 
exploiting the implications of "closed" geometries on the ground state global properties 
of the system, already studied in Ref . [3] . Such a qubit would be the elementary building 
block of a "protected" quantum computer. The task appears to be not very simple; in 
general we need a quantum system with 2^ quantum states {K being the number of big 
openings in the Josephson system under study) which are degenerate in the absence of 
external perturbations and are robust against local random fluctuations, that is against 
noise. This means that any coupling to the environment doesn't induce transitions be- 
tween the 2^ quantum states or change their relative phases. Summarizing, we need 
a system, whose Hilbert space contains a 2^-dimensional subspace characterized by the 
crucial property that any local operator O has only state-independent diagonal matrix 
elements up to vanishingly small corrections: {n \ O \m) = OoSmn + o [exp {—L)], L being 
the system size. A possible answer to such a highly non trivial requirement could be 
a system with a protected subspace built up by a topological degeneracy of the ground 
state ^15j. An alternative approach would be to exhibit a low-energy effective field theory 
for the system under study which is a topological one and whose vacua are topologically 
degenerate and, then, robust against noise. This is the approach which we follow in the 
present letter; in particular we show how to get a protected subspace with 2^ quantum 
states, = 1, by considering a Josephson junction ladder and closing it by imposing Mo- 
bius boundary conditions, in order to get a non trivial topology. We will show that such 
a system is described by a low-energy effective field theory which is a twisted conformal 
field theory [H [2] [6] [7] . Such a theory accounts very well for the topological properties of 
the system under study [3][1]. In particular we analyze the stability and transformation 
properties of the ground state wave functions under adiabatic magnetic flux change; in 
this way we are able to identify the two states of a possible protected qubit and also to 
describe its manipulation: "flip state" processes. 

We recall that Josephson junction arrays ( JJA) are a very useful tool for investigating 
quantum-mechanical behaviour in a wide range of parameters space, from Eq ^ Ej 
(where Ec = is the charging energy and Ej = is the Josephson coupling energy; 
C is the capacitance of each island and Jc is the critical current of each junction) to Ej 
^ Ec- In fact there exists a couple of conjugate quantum variables, the charge and 
phase of each superconducting island, and two dual descriptions of the array can be given 
[T7] : a) through the charges (Cooper pairs) hopping between the islands, b) through the 
vortices hopping between the plaquettes. Furthermore in the presence of an external 
magnetic field charges gain additional Aharonov-Bohm phases and, conversely, vortices 
moving around islands gain phases proportional to the average charges on the islands 
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[T8] . Such basic quantum interference effects found applications in recent proposals for 
solid state qubits for quantum computing, based on charge [19] or phase [20] degrees of 
freedom in JJAs. "Charge" devices operate in the regime Ec ^ Ej while "phase" or 
"flux" devices are characterized by strongly coupled junctions with Ej 3> Ec- 

Let us now focus on the simplest physical array one can devise in order to meet all 
the above requests, that is a Josephson junction ladder with N plaquettes closed in a 
ring geometry with a half flux quantum (|$o = threading each plaquette [5], and 
describe briefly its general properties before introducing an interaction of the charges 
(Cooper pairs) with a magnetic impurity (defect), as drawn in Fig. 1. With each site i we 
associate a phase ipi and a charge = 2enj, representing a superconducting grain coupled 
to its neighbors by Josephson couplings; rii and (pi are conjugate variables satisfying the 
usual phase-number commutation relation. The Hamiltonian describing the system is 
given by the quantum phase model (QPM): 

E / (9 \ ^ 

i V (ij) 

where Ec = -^^r- {C being the capacitance) is the charging energy at site i, while the 
second term is the Josephson coupling energy between sites i and j and the sum is over 
nearest neighbors. The most general form for the charging energy would be ^QiC^j^qj, 
where C~j'^ is the inverse capacitance matrix, but in this letter we assume for simplicity 
that the most important contribution arises from the self-energy of each grain [21] [5]. 
Aij = ^ J."' A-dl is the line integral of the vector potential associated to an external 
magnetic field B and $o = |f is the magnetic flux quantum. The gauge invariant sum 
around a plaquette is J2p^ij — ^tt/ with f = ■§^, where $ is the flux threading each 

plaquette of the ladder. Let us label the phase fields on the two legs with (p["'\ a = 1,2 
and assume Eij = E^ for horizontal links and Eij = Ey for vertical ones. Let us also make 
the gauge choice = +Tcf for the upper links, A^ = — vr/ for the lower ones and A^ = 
for the vertical ones, which corresponds to a vector potential parallel to the ladder and 
taking opposite values on upper and lower branches. 

Thus the effective quantum Hamiltonian ([T]) can be written as [S]: 



d V f d 



+ 



+ 



[2E^ cos (Xi+i - Xi) cos (0i+i - (pi- vr/) + Ey cos (20^)] , (2) 

i 

after performing the change of variables: y^-^'' = Xi + (pi, ipf'^ = Xi — (pi, where Xi, (pi 
(i.e. ^f^) are only phase deviations of each order parameter from the commensurate 
phase and should not be identified with the phases of the superconducting grains |5]. 
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Figure 1: JJL with a magnetic impurity 



When f = \ and Ec = (classical limit) the ground state of the ID frustrated 
quantum XY (FQXY) model displays - in addition to the continuous U{1) symmetry of 
the phase variables - a discrete Z2 symmetry associated with an antiferromagnetic pattern 
of plaquette chiralities Xp = il? measuring the two opposite directions of the supercurrent 
circulating in each plaquette. Thus it has two symmetric, energy degenerate, ground states 
characterized by currents circulating in the opposite directions in alternating plaquettes 
in full analogy with the checkerboard ground states of the 2D system [22]. For small Eq 
there is a gap for creation of kinks in the antiferromagnetic pattern of Xp and the ground 
state has quasi long range chiral order [5]; furthermore the charge noise, which is the 
strongest noise, has less effect in such a regime |14j|16] . 

The ladder under study can be closed and arranged in a Corbino disk geometry. As we 
will argue in the following, this is the relevant geometry for the physical implementation 
of an ideal quantum computer. In closing the ladder, we can distinguish two inequivalent 
configurations, corresponding to an even or odd number of plaquettes in the ladder. It is 
due to the antiferromagnetic pattern of plaquette chiralities, which characterizes the JJL 
ground states. 

- In the even case, the plaquettes on the opposite sides of the ladder have opposite 
chiralities, for both the degenerate ground states. So, the closed geometry can be realized 
gluing the opposite sides of the ladder keeping the ground state antiferromagnetic pattern. 

- In the odd case, the plaquettes on the opposite sides of the ladder have the same 
chiralities, for both the degenerate ground states. In this case the ladder has to be 
modified; a magnetic impurity has to be introduced, in the glue-point, which couples the 
up and down phases through its interaction with the Cooper pairs of the two legs (as 
represented in Fig.l). 

In the odd configuration, the two degenerate ground states can be represented by |0) 
and |1) and distinguished by the value of the sum over all plaquettes J2pXp respectively 
equal to —1 and +1. 
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The outlined pattern, in this closed geometry, evidences the emerging of non trivial 
topological properties intimately related to the twofold degeneracy of the ground states 
which appear to be "protected" from external perturbations [11] [1] [3]. Moreover, such 
a pattern is size independent, that is, it depends on the number of plaquettes only by 
its party and, in particular, persists also in the continuum limit oo and a — 0, 

where N is the number of plaquettes, a is the side-size of a plaquette and L = aN is the 
constant length of the ladder. Our strategy is to study the continuum limit of the J JL in 
this closed geometry; indeed, in the continuum the powerful tools of the CFT can be used 
to evince the topological properties of the system which, when extended to finite ladder, 
allow us to propose a JJL realization of a qubit device. 

In the following, we will setup the CFT analysis, while some more details on our 
twisted model (TM) can be found in the Appendix. Performing the continuum limit of 
the Hamiltonian (12]), one obtains: 



H 



dX 



+ — 



+ 



dx 



\ dx ) ~'~ ^ \dx 2 



Ey COS (20) 



(3) 



where we see that the X and fields are decoupled. In fact the X term of the above 
Hamiltonian is that of a free quantum field theory while the one coincides with the 
quantum sine-Gordon model. Through an imaginary-time path-integral formulation of 
such a model [23j it can be shown that the ID quantum problem maps into a 2D classical 
statistical mechanics system, the 2D fully frustrated XY model, where the parameter 

^ plays the role of an inverse temperature [5]. We work in the regime Ex ^ Ey 
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where the ladder is well described by our TM with central charge c = 2. 

Let us introduce in the continuum the closed geometry; in order to do so, we require 
the compactification of the ^p^""^ variables to recover the angular nature of the up and 
down fields. Then the even and odd configurations rising in the closed geometry of the 
finite ladder correspond in the continuum to two different boundary conditions for the 
fields, respectively, periodic (P) and Mobius (A) boundary conditions: 



(1) f 



0) 



X 



0) and (p\^ 



(1) 



0) 



(2) 



X 



0) 



(4) 



where we have indicated the compactified phases of the two legs as (p^^^and (p^^' , L and R 
staying for left and right components. Indeed in the limit of strong coupling the interaction 
between the magnetic impurity at point a; = (glue-point shown in Fig. 1) and the up 
and down phases gives rise to these non trivial boundary conditions for the fields [6]. 
Such a Mobius condition is naturally satisfied by the twisted field (z) of our TM (see 
eq. ( 123]) in the Appendix), which describes both the left moving component ip^^'^ and the 
right moving one (pj^ in a folded description of a system with boundary [6] [7]. In fact the 
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TM results in a chiral description of the system just described in terms of the chiral fields 
X and (see eqs. ( 1221) . (123!) ). The m-reduction technique [2] well accounts for these non 
trivial boundary conditions for the J J ladder due to the presence of a topological defect, 
already built in the construction. 

Our goal is the study of the stability and transformation properties of the four ground 
states of the JJL in the closed geometry under an adiabatic elementary flux change (±ff ) 
through the central hole of the Corbino disk. Because of the energy gap, such an adiabatic 
transformation is believed to leave the system in a ground state which can be different 
from the original one, due to the occurrence of the ground state degeneracy. This analysis 
will be the crucial step for the identification of the two states of a possible protected qubit 
and for its manipulation: "flip state" processes. 

We use the TM model to analyze such properties by standard conformal techniques. 
In the CFT description the ground state wave functions are expressed as correlation 
functions of the primary fields describing the elementary particles, in our case the Cooper 
pairs. In particular in the torus topology the characters of the theory are in one to one 
correspondence with the ground states. Indeed, as we are going to show, they describe the 
components of the "center of charge" for the corresponding ground state wave functions 
|24j . which represent coherent states of Cooper pairs on the torus. To such an extent, let 
us define for a single Cooper pair on a torus a x 6 an effective mean-field Hamiltonian 

of the kind if (x, = ifo (3;, + ^ (a;, y), where i^o (2;) y) = —ihV—2eA/c /2m is 
the Hamiltonian in the presence of an uniform magnetic field and V (x, y) is a mean-field 
scalar potential such that V {x,y) = V {x + a,y) = V {x,y + h). It is now possible to 
define the magnetic translations operators S = e^^"^^/^ and T = e^^y^l^ along the two 
cycles A and B of the torus respectively, where: 

= TT, - f % = ey = 'ny^'^BX=-%hdy^'^BX (S) 

and the gauge choice A (x, = {—By, 0) has been made. They satisfy the relations: 

S, n {x, y)] = \f, n {x, y)] =0, Sf = 62-^*^6/1-07^5 ^ (5) 

where is the magnetic flux threading the torus surface, and their action on the wave 
functions can be defined as: 

Sip [x, y) =<p{x + a,y), f<p {x, y) = e2'^*^*^V*0(^ y + h) . (7) 

Now for = M$o (i-e. when the magnetic flux is an integer number of flux quanta 



$0 = If) the condition 



5,r 



holds and we can simultaneously diagonalize the 



operators H{x,y), S, T. By introducing adimensional coordinates on the torus of the 
kind T = {uj = X + ry : X E [0, 1] , y G [0, 1]}, eqs. ([7]) can be rewritten as: 

5^ (u) = (cj + 1) , fy^ (u) = e2-^A^^(^ {00 + t) . (8) 
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One can look for eigenfunctions of the kind ^ [uS) = e™'^^ ^/ (u) and define magnetic 
translation operators Sa,Ta acting only on / (u): 

SJ {u) = f{u + a), rj {u) = e-^<"'-+2-)/ + ar) • (9) 

In this way eqs. ([H]) become: 

Sip (cu) = e'^^^y'^Sif (uj) , Tip (cu) = e'^^^y'^TJ (cu) . (10) 

Going back to the system of Cooper pairs, in order to describe a coherent state on a torus 
we look for wave functions of the kind: 

ija{u;u...,ooM) = e'''''^^"^y' fa{uJu...,UJM) , (11) 

4 

Xa Ht) , (12) 

where = is the "center of charge" variable and the non local functions Xa (^k) 

are the characters of our theory (the TM). In fact it can be shown that such characters 
are eigenfunctions of the following generalized magnetic translations operators: 

= nfll'^a/Af '^ = nfil'^V^5 (13) 

where 5*^/^/ and T^^j^ are the magnetic translation operators for the single Cooper pair. In 
this sense our characters represent highly non local functions: all the topological properties 
of our system are codified in such functions. 

On a pure topological base we expect for the torus a doubling of the ground state 
degeneracy, which can be seen at the level of the conformal blocks (characters) of our 
TM. Indeed we get for the periodic (even ladder) case an untwisted sector, P — P and 
P — A, described by the four conformal blocks (!35|) - (!38|) . and for the Mobius (odd ladder) 
case a twisted sector. A — P and A — A, described by the four conformal blocks (13T]) - (!34|) . 
Now we extract from the vacua of our theory the two states of the "protected" qubit. 

Let A be the cycle of the torus which surrounds the hole of the Corbino disk; in the 
twisted sector it is composed by the leg 1 and the leg 2 through the gluing-point as shown 
in Fig. 1. As underlined in our previous pubhcations [3j, the ground state wave functions 
of the twisted and untwisted sectors of the TM are characterized by different monodromy 
properties along the A-cycle. In particular the characters of the untwisted sector are 
single-valued functions along the A-cycle while the characters of the twisted sector pick 
up a common (—1) phase factor along the A-cycle. Such phase factors can be interpreted 
as Bohm-Aharonov phases generated while a Cooper pair is taken along the A-cycle. The 
above observation evidences a strong difference between the two inequivalent topological 
even ladder (untwisted sector) and odd ladder (twisted sector) configurations. Indeed in 
the odd ladder the ground state wave functions show a non trivial behavior implying the 
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trapping of a half flux quantum (| (|^) ) in the hole of the Corbino disk. Instead in the 
even ladder, due to the single-valued ground state wave functions, only integer numbers 
of flux quantum can be attached to the hole. 

It is worth pointing out the central role played by the "isospin" (or neutral) component 
of the TM in producing the discussed non trivial monodromy properties. To this end let us 
recall that the TM is a c = 2 CFT, composed by a c = 1 charged and a c = 1 isospin CFT 
components, as it is well evidenced by the character decompositions in the Appendix. 
Now the transport of a Cooper pair along the A-cycle is implemented by a simultaneous 
and identical translation = Aw„ = 2 of the charged and the isospin variables. The 
"charged" characters have trivial monodromy with respect to this transformation, being: 

Ki{w, + 2\t)=Ki{w,\t), / = 0,..,3, (14) 

while the "isospin" contribution is the one responsible for the non trivial monodromy of 
the complete ground state wave functions: 

Xo,i(2|r) = Xo,i(0|r) , X^S'^\r) = (-l)x^(0|r) (15) 

and the same is true for the characters X/s- Let us notice that the change in sign in 
the last relation of eq. f|T5|) shows the presence in the spectrum of excitations carrying 
fractionalized charge quanta. More precisely the presence in the isospin component of one 
twist-field (with conformal dimension A = 1/16) characterizes all the conformal blocks 
of the twisted sector and accounts for the trapping of a half flux quantum in the hole 
of the closed JJL. At this point it is crucial to observe that in order to create such a 
fractionally charged excitation in the ground state a finite energy must be provided, so 
assuring the presence of a finite gap separating the ground state from the excited ones 
(that is in complete analogy with the presence of a gap separating the ground state from 
higher energy states in the Laughlin Hall fluid pS]). 

We are now in the position to address the study of the stability and transformation 
properties of the ground state wave functions when a magnetic flux change takes place 
through the central hole of the closed JJL. The above analysis shows that at the level 
of the wave functions it has the effect to change the monodromy along the A-cycle due 
to the corresponding change in the Bohm-Aharonov phase. Such a modification can be 
implemented on the center of charge component of the wave function, i.e. the characters, 
with a well defined transformation. In the case of the charged component this analysis 
has been brought out already in in the physical contest of the quantum Hall effect. 
Let us adapt here the results for the charged component of our TM. 

On a pure physical ground the fact that we are considering a magnetic flux change, 
which is on one side integer in the flux quantum (one flux quantum change ±^) and on 
the other side adiabatic, suggests both that the monodromy properties do not change and 
that the system remains in a degenerate ground state. Such a physical picture is in fact 
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confirmed for the charged component of our TM; indeed, the flux change is implemented 
on the level of charged characters by the transformation ^/2- 

^,Ki{w,\T)=e("^y'''^^"-^Ki(w, + ^\r'^ = i^m K|r) , / = 0,..,3. (16) 

In particular the charged component wave functions realize a flip process {I —>■ l + l) under 
one magnetic flux quantum change. 

However the analysis for the complete TM, with charged and isospin components, is 
more involved. In particular the problem of the stability of the ground state wave func- 
tions under the change of one flux quantum in the central hole has to be clearly brought 
out. This is mainly due to the non trivial interplay between charged and isospin com- 
ponents summarized in the so-called m-ality parity rule, which characterizes the gluing 
condition for the charged and isospin excitations (see Appendix). The main point being 
the compatibility between such parity rule and the transformation of the complete char- 
acters of TM under the insertion of a flux quantum through the hole of the closed ladder, 
which reads as: 

Ti/2f{wMr) = e2-("'-+-(-"+-=))/(^„ + aT\w, + ar|r) , (17) 

a=l/2 

where f{wn = 0\wc\t) stays for any character of TM. The full list of such transformations 
is presented in Appendix, here we only comment on the very simple and clear picture 
which emerges for the stability and transformations of the ground states of the closed 
JJL. 

The even configuration of the closed JJL (periodic case) is proven to be unstable under 
this transformation. Indeed eq. fHOl) and (j^Tl) show the decoupling of the untwisted P — A 
sector and of the state Xa of the P — P sector while eq. P2l) shows that the state Xp of 
P — P sector gets excited by this transformation. 

For the odd configuration of the closed JJL (Mobius case), eq. fH3|) shows that the 
twisted A — A sector decouples. So we are left only with the A — P sector, with the two 
ground states flipping one into the other under an adiabatic flux change of through 
the central hole, as it can be seen from eq. (Hll). 

Summarizing, between the two inequivalent configurations for the closed JJL, corre- 
sponding to even and odd ladder, just the odd one is proven to be stable under an adiabatic 
flux change of through the central hole. Then in terms of the ground states center 
of charge wave functions (characters), we can make the following identiflcations: 

|0)-xJ)(OK|r), |l)^x^)(OK|r). (18) 

Then |0) and |1) are the two ground states of the odd closed JJL characterized by the 
size invariant sum over all plaquettes YlpXp respectively equal to —1 and -|-1. 
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(a) (b) 

Figure 2: The two logical states of the JJL qubit: a) for the state |1), b) for the state |0) . 

Based on the above consideration we are ready to propose the odd closed JJL as our 
protected qubit. In fact the two ground states |0) and |1) work as the two logical states 
of the qubit and the required one qubit operations: 

10)-. |1), |1)-|0), 

are simply implemented by insertion of a flux quantum (±|f ) through the central hole. 

Let us now make a comment on the stability of such qubit device under local pertur- 
bations. Local perturbations can be viewed in such a context as creating a finite energy 
excitation above the ground states in the form of double kinks. A double kink can be 
produced from the ground state by exchanging the chirality of two nearest neighbor pla- 
quettes in the ladder and, as such, it is local and it leaves invariant the chirality sum 

Xp over all plaquettes and so the characterization of the two logical states. Further- 
more, since a double kink can be described by the presence of two elementary half flux 
quanta of opposite sign (±Hf ) localized in between the pairs of plaquettes with the same 
chirality, it doesn't produce any flux change in the central hole. In this way the excited 
logical state wave function shows the same monodromy properties along the A-cycle as 
the corresponding ground state one and, in particular, satisfies the same transformation 
rules under an adiabatic flux change of through the central hole. Summarizing, the 
characterization of the two logical states and their flipping processes are left unchanged 
under local perturbations, which produce a finite energy excitation above the ground 
state. 

The minimal configuration for such a protected qubit is represented in Fig. 2 by a 
closed fully frustrated JJL with = 3 plaquettes, 3 being the minimum odd number of 
plaquettes needed in order to fulfill all the above requests. 

Now it should be possible to construct symmetric (s) and antisymmetric (a) linear 
combinations of such degenerate ground states and then to control their amplitude and 
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relative phase: such operations are needed in order to prepare the qubit in a definite 
state and to manipulate it [27] . In order to realize the logical NOT we must perform an 
adiabatic change of local magnetic fields that drags one vortex across the system, i.e. a 
fiux quantum through the A-cycle of the torus, and flips the state of the system, so lifting 
the degeneracy [Hj . 

Josephson junction ladders with annular geometry have been fabricated within the 
trilayer Nh/Al — AlOx/Nb technology and experimentally investigated [28,J, but in such 
a case the application of an external transverse magnetic field is needed in order to fulfill 
the requirement of full frustration and that could be another source of decoherence. It is 
now possible to avoid such a problem by realizing arrays with a built-in frustration. In 
fact high-Tc Josephson junction arrays have been recently proposed [22], which support 
degenerate spontaneous current states in zero magnetic field due to the presence of pla- 
quettes containing an odd number of vr-junctions [30]. Such unconventional junctions can 
be realized because of the d-wave symmetry of high-Tc superconductors [31] , which pro- 
duces a TT-shift in the phase of the wave function on one side of the junction. Furthermore 
TT-junctions can be obtained also with superconducting-ferromagnetic-superconducting 
junctions (SFS) [52] • In this way it is possible to avoid the external frustration bias but, 
in any case, external magnetic fields are needed for control and read-out operations: in 
fact our JJL qubit is a fiux device [20] . 

So in principle an experimental setup for the realization of our protected qubit can be 
conceived: the JJL, arranged in a ring geometry (Corbino disk) with an odd number of 
plaquettes along the inner hole, should be equipped with a coil Phias which can be used 
to set the system in one of the two ground states |0), |1); another coil Phoie is needed 
in order to control the mixing of the two ground states and to carry out the fiipping. 
Finally read-out coils, coupled to external conventional SQUIDs, are needed in order 
to read out the state of the system after the quantum evolution. The whole device has 
to be embedded in a superconducting cavity in order to guarantee the stability of the 
boundary conditions. In this way a reliable qubit is built up, whose quantum evolution 
can be controlled in order to perform all possible single qubit logical operations. Then 
such qubits can be antiferromagnetically coupled by means of suitable superconducting 
fiux transformers ( "eight" coils) which provide an inductive coupling and whose strength 
can be controlled within a wide range of useful values: in this way a quantum register can 
be realized and all multi-qubits logical operations can be performed. 

Summarizing, a single non-interacting qubit is described by a double well potential 
and the external magnetic fiux controls the energy difference between the minima, the 
symmetric situation being for = 0. Each logical state, |0) or |1), is represented by a 
wave function localized in a distinct potential well and corresponds to distinguishable fiux 
states trapped in the plaquettes of the ladder with current fiowing in opposite directions 
in alternating plaquettes. When the energy difference e of the minima of the two different 
wells is small with respect to the oscillation frequency u around the minima, e <^ cu, these 
states become coupled and the wave functions spread over both the wells, the coupling 
being maximum in resonance conditions (e = 0), while the energy eigenstates tend to be 
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localized in one of the well away from resonance as e is increasing. The coupling of the 
states can be described by a tunneling amplitude A (e) and the effective Hamiltonian of 
any qubit reduces to the regular two-state form in the basis of these logical states: 



Heff = \[e (|0) (0| - |1) (1|) - A (|0) (1| + |1) (0|)] = \ {ea^ - Aa^ , (19) 



where a^, are the Pauli spin matrices. 

The diagonal elements of H^ff can be easily controlled by an external magnetic field in 
the z-direction producing an external flux while the off-diagonal elements are related to 
the tunneling amplitude and thus are controlled by the adiabatic change of the magnetic 
flux in the central hole. The general state vector of such a qubit is the linear combination 
of the basis states: 



so it is described by two complex numbers a and (3. 

When the inductive coupling among qubits is turned on, there could be a bias in, 
say, the j-th qubit even though = and, as a consequence, its logical states may be 
asymmetric. In the approximation in which every JJL can be considered as a two level 
system, the system of flux linked qubits can be described by an effective Hamiltonian of 
the kind: 



In order to control such Hamiltonian, one should be able to modulate the tunneling 
amplitude of each qubit as well as to switch on and off the magnetic coupling between 
neighbors qubits. The analysis of multi-qubit logical operations will be the subject of a 
future publication. 

In conclusion in this letter we have presented a simple collective description of a 
fully frustrated ladder of Josephson junctions arranged in a non trivial geometry, with 
a macroscopic half flux quantum trapped in the hole. The powerful tools of the CFT 
have been used to evince the topology, the stability and the transformation properties of 
the system. In particular it has been shown how such features can be exploited for the 
realization of a "protected" qubit: a simple device has been proposed and its operation 
mode has been briefly sketched. 
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Appendix: reminder of TM for the JJL 



Here we briefly summarize the main results of our theory, the TM, for the fully frustrated 
JJL [1][3]. We first construct the bosonic theory and show that its energy momentum 
tensor fully reproduces the Hamiltonian of eq. Q for the JJL. That allows us to describe 
the JJL excitations in terms of the primary fields Va (z). Then we show that it is possible 
to construct the A^— vertices correlator for the torus topology in 2D (basically by letting 
the edge to evolve in "time" and to interact with external vertex operators placed at 
different points). We assume that a suitable correlator is apt to describe the ground 
state wave function of the JJL at T = temperature and then perform an analysis of 
the symmetry properties of its center of charge wave function (conformal blocks), which 
emerge in the presence of vortices carrying half quantum of flux (ff))- 

Let us focus on the m-reduction procedure |2j for the special m = 2 case (see Ref. [1] 
for the general case), since we are interested in a system with a Z2 symmetry and choose 
the "bosonic" theory [1][3], which well adapts to the description of a system with Cooper 
pairs of electric charge 2e in the presence of a topological defect [6], i.e. a fully frustrated 
JJL. To each of the two legs (edges) of the ladder we assign a chirality, so making a 
correspondence between up-down leg and left-right chirality states. 

Let us now write each phase fleld as the sum ip^°-'> {x) = y^^^ (x) + Lp^^ (x) of left and 
right moving flelds deflned on the half-line because of the topological defect located in 
X = 0. Then let us deflne for each leg the two chiral flelds Lp^e)> {x) = Lp^l^ (x) ± Lp'^^ {—x), 
each deflned on the whole x— axis [35]. In such a framework the dual flelds Lp'o'^ (x) are 
fully decoupled because the corresponding boundary interaction term in the Hamiltonian 
does not involve them [36]; they are involved in the deflnition of the conjugate momenta 

H(a) = {^x'^'o'^^ = (^^"T^) present in the quantum Hamiltonian. Performing the change 

of variables if^e^ = X + (p, (pf^ = X — (p {(f^o^ = X + (p, Lp'if' = X — for the dual ones) 
we get the quantum Hamiltonian ([3]) but now all the flelds are chiral ones. Finally let 
us identify in the continuum such chiral phase flelds a = 1,2, each deflned on the 
corresponding leg, with the two chiral flelds Q^°'\ a = 1, 2 of the TM with central charge 
c = 2. 

As a result of the 2-reduction procedure [1] [2j we get a c = 2 orbifold CFT, the TM, 
whose flelds have well deflned transformation properties under the discrete Z2 (twist) 
group, which is a symmetry of the TM. Its primary flelds content can be expressed in 
terms of a Z2-invariant scalar fleld X{z), given by 

X(^) = i(g«(^)+Q(2)(^))^ (22) 
describing the continuous phase sector of the theory, and a twisted fleld 

<P{z) = l{Q^'\z)-Q^'Hz)), (23) 
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which satisfies the twisted boundary conditions (j){e^'^z) = —4>{z) More exphcitly such 
a field can be written in terms of the left and right moving components ^p^'^ as 
we stated above; then the Mobius boundary conditions given in eq. (jl]) are described 
by the boundary conditions for 0. This will be more evident for closed geometries, i.e. 
for the torus case, where the magnetic impurity gives rise to a line defect in the bulk, 
so allowing us to resort to the folding procedure and introduce boundary states [6] [7]. 
Such a procedure is used in the literature to map a problem with a defect line (as a bulk 
property) into a boundary one, where the defect line appears as a boundary state of a 
theory which is not anymore chiral and its fields are defined in a reduced region which is 
one half of the original one. Our approach, the TM, is a chiral description of that, where 
the chiral field defined in (— L/2, L/2) describes both the left moving component and 
the right moving one defined in (—L/2, 0), (0, L/2) respectively, in the folded description 
[B][7]. Furthermore to make a connection with the TM we consider more general gluing 
conditions: 



the — (+) sign staying for the twisted (untwisted) sector. We are then allowed to use the 
boundary states given in [37] for the c = 1 orbifold at the Ising^ radius. The X field, which 
is even under the folding procedure, does not suffer any change in boundary conditions [S] 
while condition (jl]) is naturally satisfied by the twisted field {z). So topological order can 
be discussed referring to the characters with the implicit relation to the different boundary 
states (BS) present in the system [6]. These BS should be associated to different kinds of 
linear defects compatible with conformal invariance. 

The fields in eqs. f l22p - (l25]) coincide with the ones introduced in eq. ([2]). In fact the 
energy momentum tensor for such fields fully reproduces the second quantized Hamilto- 
nian of eq. ([3]). The whole TM theory decomposes into a tensor product of two CFTs, 
a twisted invariant one with c = | and the remaining c = | one realized by a Major ana 
fermion in the twisted sector. In the c = | sub-theory the primary fields are composite 
vertex operators V {z) = U'^ {z) ip {z) or Vqh {z) = U'^ [z) a {z), where 



is the vertex of the continuous sector with ai = |, / = 1,...,4 for the SU{2) Cooper 
pairing symmetry used here. Regarding the other component, the highest weight state in 
the isospin sector, it can be classified by the two chiral operators: 



twist) or Neveu-Schwartz {Z2 twisted) boundary conditions [T][2] in a fermionized version 
of the theory. The Ramond fields are the degrees of freedom which survive after the 
tunnelling and the parity symmetry, which exchanges the two Ising fermions, is broken. 



0l(x = 0) = t<Pr{x = 0) 
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Besides the fields appearing in eq. ( l25ll . there are the <J (z) fields, also called the twist 
fields, which appear in the quasi-hole primary fields Vqh (z) . The twist fields have non 
local properties and decide also for the non trivial properties of the vacuum state, which 
in fact can be twisted or not in our formalism. 

Starting from the primary fields {z) we can now construct the non perturbative 
ground state wave function of the JJL system for the torus topology. It turns out that 
by construction it results as a coherent superposition of gaussian states with all the non 
trivial global properties of the order parameter. In fact by using standard conformal field 
theory techniques it is now possible to generate the torus topology, starting from the edge 
theory, just defined above. That is realized by evaluating the A^-vertices correlator 

(n|V;(;2i)...V;(;2^)e2--^°|n), (26) 

where Va (zi) is the generic primary field representing the excitation at Zi, Lq is the 
Virasoro generator for dilatations and r the proper time. The neutrality condition X] = 
must be satisfied and the sum over the complete set of states \n) is indicating that a trace 
must be taken. It is very illuminating for the non expert reader to pictorially represent 
the above operation in terms of an edge state (that is a primary state defined at a given 
r) which propagates interacting with external fields aX zi . . . zn and finally getting back to 
itself. In such a way a 2D surface is generated with the torus topology. It is interesting to 
observe that such a procedure is equivalent to the coherent insertion of correlated relevant 
vortices (as provided by the CFT description) at positions Zi . . . Z]\f , as they appear in 
the non perturbative ground state of the physical JJL system. From such a picture it is 
evident then how the degeneracy of the non perturbative ground state is closely related 
to the number of primary states. Furthermore, since in this letter we are interested in the 
understanding of the topological properties of the system, we can consider only the center 
of charge contribution in the above correlator, so neglecting its short distances properties. 
To such an extent the one-point functions are extensively reported in the following. 

On the torus [2J the TM primary fields are described in terms of the conformal blocks 
of the Z2-invariant c = | sub-theory and of the non invariant c = | Ising model, so 
refiecting the decomposition on the plane above outlined. The characters Xo(0|t), Xi(0|T), 
XJ-(0|t) express the primary fields content of the Ising model [M] with Neveu-Schwartz 
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{Z2 twisted) boundary conditions [2], while 

X(ot^'(0|t/7,|r) = xo{0\T)Koiw,\T) + xi^{0\T)K2{w,\T), (27) 

xI;'^\0\w,\t) = x±{0\r){K,{w,\r) + K,{w,\T)), (28) 

x";'^\0\w,\t) = xi{0\T)Ko{w,\r) + xo{0\T)K2{w,\T) (29) 

represent those of the Z2-invariant c = | CFT. They are given in terms of a "charged" 
Kaiwrlr) contribution: 



K2i+i{w\T) = — ^ e 

r/(r) 



2l+i 
4 



(2u;|4r) , with / = 0, 1 and i = 0, 1 , (30) 
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and a "isospin" one X/3(0|t), (the conformal blocks of the Ising Model), where Wc = 

In Zc is the torus variable of the "charged" component while the corresponding argu- 

27rz 

ment of the isospin block is w„ = everywhere. 

If we now turn to the whole c = 2 theory, the characters of the twisted sector are given 

by: 

xJ)(OK|r) = x^(0|r)(xo + Xi)(0|r)(Ko + i^2)(«^c|r), (31) 
X+)(OK|r) = x^(0|r)(xo + Xi)(0|r)(Ki + ir3)(^c|r), (32) 

for the A — P sector and by: 

X^o)(Okc|r) = x±^{0\T)(xo-Xi){0\r){Ko-K2){w,\T), (33) 
X(i)(0|«;e|r) = x±.{0\r)(xo-Xi){0\r)iK, + K,){w,\T), (34) 

for the A — A one. Furthermore the characters of the untwisted sector are [2]: 

Xlo)iO\wc\T) = (^xoXo - XiXi) (0|r)i^o(wc|r) + (xoXi - XiXo) (0|r)i^2 (^dr), (35) 

X^i)(0|wc|r) = (xoXi - XiXo) {0\r)Koiwc\r) + (^xoXo - XiXi) (0|r)ir2 (wc|r), (36) 
for the P — A sector while for the P — P sector we have: 

xMwc\r) = l{xo- Xi) (0|r) (xo - Xi) (0|r)(iro - K2){w,\t) , (37) 

X+(0|we|r) = l{xo + Xi) (0|r) (xo + Xi) (0|r)(iro + i^2)(«^c|r), (38) 



and 



X+(0|t/7e|r) = X^(0|r)x 1 (0|r) (K, + K,) («;,|r). (39) 



Let us comment that the above factorization expresses the parity selection rule (m-ality), 
which gives a gluing condition for the "charged" and "isospin" excitations. 

It is worth underlining that in the P — P sector, unlike for the other sectors, modular 
invariance constraint requires the presence of three different characters. The isospin 
operator content of the character x+(0|wc|r) clearly evidences its peculiarity with respect 
to the other states of the periodic (even ladder) case. Indeed it is characterized by two 
twist fields (A = 1/16) in the isospin components. The occurrence of the double twist 
in the state described by (0|wc|t) is simply the reason why such a state is a periodic 
state. Indeed, being an isospin twist field the representation in the continuum limit of a 
magnetic impurity (a half flux quantum trapping or equivalently a kink), the double twist 
corresponds to a double half flux quantum trapping, i.e. one flux quantum, typical of the 
periodic configuration. 
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The above analysis would suggest that the P—P state described by (0|uvV_;|r) embeds 
in the continuum limit a kink-antikink excitation, i.e. it represents an excited state in the 
P — P sector. In this way, as it happens for all the other sectors, the P — P sector is left 
with just two degenerate ground states ( x+(0|wc|t) and x^(0|wc|t)) and, as expected on 
a pure topological base, the ground state degeneracy in the torus topology is the double 
of that of the disk. 

Let us now present the full list of character transformations under the insertion of a 
magnetic flux quantum through the hole of the closed ladder. 

In the even closed JJ ladder configuration, we have that the two ground state wave 
functions of the P — A sector decouple, being 

Ti/2X^O)(0|We|r)=0, Ty2Xi,){0\Wc\T) = 0. (40) 

Concerning the P — P sector, we have: 

ri/2x+(okc|T) = (41) 

and 

Ti/2xJ(0|we|r) = x+(Okc|r) ( Ti/2X+(0|w;c|r) = X^Olwdr) ). (42) 

Such transformations show the instability of the P — P sector under the insertion of a flux 
quantum through the hole of the closed ladder. More precisely the state Xa (0|'?^c|t) decou- 
ples while the state x^(0|wc|t) gets excited to the state with a kink-antikink configuration 
X+(0|wc|t). 

Furthermore in the odd closed JJ ladder configuration, we have that the two ground 
state wave functions of the A — A sector decouple, being 

^i/2X^o)(0kc|r) = 0, Ti/2X^,)(0|w,|r) = 0. (43) 

Concerning the A — P sector, we have that the two ground state wave functions transform 
as: 

Ti/2x5)(0kc|r) = xfi)(Oke|r) , ri/2x5)(0|w,|r) = xto)i^\wc\r) . (44) 

Concluding, the full set of transformations, here presented, allows to claim the fol- 
lowing simple and clear picture: the odd closed JJL configuration is the only one which 
is stable under the insertion of a magnetic flux quantum through the central hole; more- 
over, in such odd JJL configuration such a magnetic fl,ux insertion simply implements the 
flipping process between the two degenerate ground states |0) and |1). 

References 

[1] G. Cristofano, G. Maiella, V. Marotta, Mod. Phys. Lett. A 15 (2000) 1679. 



18 



[2] G. Cristofano, G. Maiella, V. Marotta, G. Niccoli, Nucl. Phys. B 641 (2002) 547. 

[3] G. Cristofano, V. Marotta, A. Naddeo, G. Niccoli, Eur. Phys. J. B 49 (2006) 83. 

[4] G. Cristofano, V. Marotta, A. Naddeo, J. Stat. Mech.: Theor. Exper. (2005) P03006. 

[5] M. Kardar, Phys. Rev. B 30 (1984) 6368; M. Kardar, Phys. Rev. B 33 (1986) 3125; 
E. Granato, Phys. Rev. B 42 (1990) 4797: E. Granato, J. Appl. Phys. 75 (1994) 
6960; E. Granato, Phys. Rev. B 72 (2005) 104521. 

[6] G. Cristofano, V. Marotta, A. Naddeo, Phys. Lett. B 571 (2003) 250. 

[7] G. Cristofano, V. Marotta, A. Naddeo, Nucl. Phys. B 679 (2004) 621. 

[8] G. Cristofano, V. Marotta, A. Naddeo, G. Niccoli, J. Stat. Mech.: Theor. Exper. 
(2006) L05002. 

[9] O. Foda, Nucl. Phys. B 300 (1988) 611. 

[10] G. Cristofano, V. Marotta, P. Minnhagen, A. Naddeo, G. Niccoli, J. Stat. Mech.: 
Theor. Exper. (2006) P11009. 

[11] X. G. Wen, Q. Niu, Phys. Rev. B 41 (1990) 9377. 

[12] X. G. Wen, Phys. Lett. A 300 (2002) 175; X. G. Wen, Phys. Rev. B 65 (2002) 
165113; T. Senthil, M. P. A. Fisher, Phys. Rev. B 62 (2000) 7850; T. Senthil, M. P. 
A. Fisher, Phys. Rev. B 63 (2001) 134521. 

[13] M. C. Diamantini, P. Sodano, C. A. Trugenberger, Nucl. Phys. B 474 (1996) 641; 
M. C. Diamantini, P. Sodano, C. A. Trugenberger, J. Phys. A 39 (2006) 253; M. C. 
Diamantini, P. Sodano, C. A. Trugenberger, Eur. Phys. J. B 53 (2006) 19. 

[14] L. B. loffe, V. B. Geshkenbein, M. V. Feigelman, A. L. Fauchere, G. Blatter, Nature 
398 (1999) 679; G. Blatter, V. B. Geshkenbein, L. B. loffc, Phys. Rev. B 63 (2001) 
174511; L. B. loffe, M. V. Feigelman, Phys. Rev. B 66 (2002) 224503; B. Doucot, J. 
Vidal, Phys. Rev. Lett. 88 (2002) 227005; B. Doucot, M. V. Feigelman, L. B. loffe, 
Phys. Rev. Lett. 90 (2003) 107003; L. B. loffe, M. V. Feigelman, A. S. loselevich, D. 
Ivanov, M. Troyer, G. Blatter, Nature 415 (2002) 503. 

[15] A. Kitaev, Ann. Phys. 303 (2003) 2. 

[16] B. Doucot, L. B. loffe, J. Vidal, Phys. Rev. B 69 (2004) 214501; B. Doucot, M. V. 
Feigelman, L. B. loffe, A. S. loselevich, Phys. Rev. B 71 (2005) 024505; B. Doucot, 
L. B. loffe, Phys. Rev. A 72 (2005) 032303; B. Doucot, L. B. loffe. New J. Phys. 7 
(2005) 187; B. Doucot, L. B. loffe, 0704.0900vl. 

[17] Y. M. Blanter, R. Fazio, G. Schon, Nucl. Phys. B S58 (1997) 79. 



19 



[18] Y. Aharonov, A. Casher, Phys. Rev. Lett 53 (1984) 319. 

[19] A. Shnirman, G. Schon, Z. Hermon, Phys. Rev. Lett. 79 (1997) 2371. 

[20] T. P. Orlando, J. E. Mooij, L. Tian, C. H. van der Wal,L. S. Levitov, S. Lloyd, J. J. 
Mazo, Phys. Rev. B 60 (1999) 15398. 

[21] R. M. Bradley, S. Doniach, Phys. Rev. B 30 (1984) 1138. 

[22] S. Teitel, C. Jayaprakash, Phys. Rev. Lett. 51 (1983) 1999; S. Teitel, C. Jayaprakash, 
Phys. Rev. B 27 (1983) 598. 

[23] J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, Oxford Science Pub- 
lications, Oxford (1989). 

[24] G. Cristofano, G. Maiella, R. Musto, F. Nicodemi, Nucl. Phys. B 33C Proc. Suppl. 
(1993) 119. 

[25] R. B. Laughlin, Phys. Rev. B 23 (1981) 5632; R. B. Laughlin, Phys. Rev. Lett. 50 
(1983) 1395; R. B. Laughlin, Rev. Mod. Phys. 71 (1999) 863. 

[26] G. Cristofano, G. Maiella, R. Musto, N. Nicodemi, Mod. Phys. Lett. A 7 (1992) 
2583. 

[27] A. Ekert, R. Jozsa, Rev. Mod. Phys. 68 (1996) 733; A. Steane, Rep. Prog. Phys. 61 
(1998) 117; A. Galindo, M. A. Martin-Delgado, Rev. Mod. Phys. 74 (2002) 347. 

[28] P. Binder, D. Abraimov, A. V. Ustinov, S. Flach, Y. Zolotaryuk, Phys. Rev. Lett. 84 
(2000) 745; P. Binder, A. V. Ustinov, Phys. Rev. E 66 (2002) 016603. 

[29] H. Hilgenkamp, H. J. Ariando, H. Smilde, D. H. A. Blank, G. Rijnders, H. Rogalla, 
J. R. Kirtley, C. C. Tsuei, Nature 422 (2003) 50; J. R. Kirtley, C. C. Tsuei, H. J. 
Ariando, H. Smilde, H. Hilgenkamp, Phys. Rev. B 72 (2005) 214521. 

[30] A. I. Larkin, Yu. N. Ovchinnikov, Sov. Phys. Rev. JETP 20 (1965) 762; P. Fulde, R. 
A. Ferrel, Phys. Rev. 135 (1964) A550. 

[31] J. R. Kirtley, C. C. Tsuei, Rev. Mod. Phys. 72 (2002) 969; A. A. Golubov, M. Yu. 
Kuprianov, E. Il'ichev, Rev. Mod. Phys. 76 (2004) . 

[32] V. V. Ryazanov, V. A. Oboznov, A. Y. Rusanov, A. V. Veretennikov, A. A. Golubov, 
J. Aarts, Phys. Rev. Lett. 86 (2001) 2427; V. V. Ryazanov, V. A. Oboznov, A. V. 
Veretennikov, A. Y. Rusanov, Phys. Rev. B 65 (2002) 020501. 

[33] R. K. BuUough, N. M. Bogoliubov, G. D. Pang, J. Timonen, in Chaos, Solitons and 
Fractals, ed. M. Lakshmanan, 5 (1995) 2639. 

[34] N. M. Bogoliubov, A. G. Izergin, N. A. Kitanine, Nucl. Phys. B 516 (1998) 501. 



20 



[35] F. D. M. Haldanc, Phys. Rev. Lett. 47 (1981) 1840; E. Orignac, T. Giamarchi, Phys. 
Rev. B 57 (1998) 11713; E. Orignac, T. Giamarchi, Phys. Rev. B 64 (2001) 144515. 

[36] I. Affleck, A.W.W. Ludwig, Nucl. Phys. B 352 (1991) 849; I. Affleck, A.W.W. Lud- 
wig, Nucl. Phys. B 360 (1991) 641; I. Affleck, Acta Phys. Pol. 26 (1995) 1869. 

[37] M. Oshikawa and I. Affleck, Nucl. Phys. B 495 (1997) 533. 

[38] P. Di Francesco, P. Mathieu and D. Senechal, Conformal Field Theories, Springer- 
Verlag, Berlin (1996). 



21 



